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Abstract 

There exist certain intrinsic relations between the ultraviolet divergent graphs 
and the convergent ones at the same loop order in renormalizable quantum 
field theories. Whereupon we may establish a new method, the intrinsic 
regularization method, to regularize those divergent graphs. In this paper, 
we apply this method to QCD at the one loop order. It turns out to be 
satisfactory: The gauge invariance is preserved manifestly and the results 
are the same as those derived by means of other regularization methods. 
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1 Introduction 



Over the decades, as is well known, a wide variety of regularization schemes have been developed 
in quantum field theory [H]. However, as every schemes heve their own distinct advantages and 
disadvantages, this topic is still one of the important and fundamental issues under investigation. 
One of the most challenging problem is perhaps how to preserve all properties of the original action 
manifestly and consistently. 

A few years ago, a new regularization method named intrinsic vertex regularization was first 
proposed for the ^ 4 theory by Wang and Guo [||. The key point of the method is, in fact, based 
upon the following simple observation: For a given ultraviolet divergent function at certain loop 
order in a renormalizable QFT, there always exists a set of convergent functions at the same loop 
order such that their Feynman graphs share the same loop skeleton and the main difference is 
that the convergent ones have additional vertices of certain kind and the original one is the case 
without these vertices. This is, in fact, a certain intrinsic relation between the original ultraviolet 
divergent graph and the convergent ones in the QFT. It is this relation that indicates it is possible 
to introduce the regularized function for the divergent function with the help of those convergent 
ones so that the potentially divergent integral of the graph can be rendered finite while for the 
limiting case of the number of the additional vertices q — > the divergence again becomes manifest 
in pole(s) of q. 

To be concrete, let us consider a 1PI graph with / internal lines at one loop order in the 4 
theory. Its superficial degree of divergences in the momentum space is 

5 = 4-21. 

When I = 1 or 2, the graph is divergent. Obviously, there exists such kind of graphs that they 
have additional q four-i^-vertices in the internal lines. Then the number of internal lines in these 
graphs is I + q so that the divergent degree of the new 1PI graphs become 

5' = 4-2(1 + q). 

If q is large enough, the new ones are convergent and the original divergent one is the case of q = 0. 
Thus, a certain intrinsic relation has been reached between the original divergent 1PI graph and 
the new convergent ones at the same loop order. 

However, application of this method to QED runs into a difficulty. The problem is that, unlike 
the the cj> theory, the electron-photon vertex in QED carries a 7-matrix and is a Lorentz vector. 
As a result, simply inserting the vertex would increase the rank of the function as Lorentz tensors 
and would make the problem quite complicated. In order to overcome this difficulty, in §, [§ the 
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authors introduced an alternative method. We follow the example of the </> 4 theory to demonstrate 
the key point of this method: One shifts the mass term of the cj) field from m 2 to m 2 +fj? and regards 
[x 2 cp 2 as a new vertex in addition to the vertex Xcj> . By inserting the new vertex into the internal cj> 
lines in the graph of a given 1PI n-point divergent function, a set of new convergent functions can 
be obtained provided the number of inserted vertices, q, is large enough. Then one can introduce 
a new convergent function, the regularized function, and the potential infinity in the original 1PI 
n-point function may be recovered as the q — > limiting case of that function. Obviously, the 
mass shifting method can be easily generalized to QED by simply shifting the electron mass from 
mtom||i and regarding the term —fiipip as a new vertex. In fact, as has been shown in Q, it 
turns out to be successful to QED. Nevertheless, it is not really intrinsic since the Feynman rules 
of the theory have to be modified. As a result, it may not completely work for non-Abelian gauge 
theories, e.g., QCD, because generally it is not clear whether the gauge symmetry can be preserved 
manifestly for these theories, although such a proof for QED at one loop level has been given ||. 

Very recently, we presented an improved approach in |7| to reexamine the cj) 4 theory and QED, 
in which a new concept, inserter, was introduced. An inserter is a vertex or a pair of vertices linked 
by an internal line, in which the momenta of the external legs are all set to zero, and, if there are 
any, all the Lorentz indices are contracted in pair by the spacetime metric and all the internal gauge 
symmetry indices are contracted by the Killing-Cartan metric in the corresponding representation, 
so that as a whole an inserter always carries the vacuum quantum numbers, i.e. zero momentum, 
scalar in the spacetime symmetry, and singlet in internal and gauge symmetries. It is not hard to 
see that in any given QFT as long as a suitable kind of inserters are constructed with the help 
of the Feynman rules of the theory, some intrinsic relations between the divergent functions and 
convergent ones at the same loop order will be found by simply regarding the convergent ones as 
the ones given by suitably 

inserting q-inserters in all internal lines in the given divergent ones. The crucial point of this 
approach, therefore, is very simple but fundamental, that is, the entire procedure is intrinsic in the 
QFT. There is nothing changed, the action, the Feynman rules, the spacetime dimensions etc. are 
all the same as that in the given QFT. This is a very important property which should shed light 
on the challenging problem mentioned at the beginning of the paper. Consequently, in applying to 
other cases all symmetries and topological properties there should be preserved in principle. 

In what follows, we concentrate on how to apply the inserter approach to QCD at one loop 
order. We present the main steps and the results of the inserter regularization procedure for it. 
We find that, as is expected, the gauge invar iance is preserved manifestly, and all results are the 
same as those derived by means of other regularization methods. 



3 



2 Intrinsic Regularization in QCD 



The QCD Lagrangian, including ghost fields and gauge fixing terms, can be written as 

Cqcd = -\{d»A a v - d v Al + g c f abc AlA c v ){d»A av - d u A a ^ + g c f ade A d »A e ») 

-^(d^) 2 - td^5 ac - g c f abc A b »)r, c + #7^ " fcc^y) - mty , (1) 

and the Feynman rules are well known. 

The main steps of the inserter approach for QCD may be stated more concretely as follows. 
First, we should construct the inserters in QCD. This work, with regards to simplicity and consis- 
tency with other theories, e.g., the electroweak theory, may actually be done within a more general 
framework, namely, within the framework of the standard model in which QCD is contained. The 
explicit expressions of all inserters in the standard model have been preestablished in Jjj. Here we 
merely list those relevant to QCD: 



The gluon- inserter: 



I {9 i?(p) = Sig 2 c C 2 (8)g, u 5 ab . 



The ghost-inserter: 

I {9f iaM = -ig 2 c C 2 (8)5 aia2 . (3) 



• The quark-inserter: 

I {Q} (P) = (4) 

In eqs.@ and (||), 6*2(8) is the second Casimir operator valued in the adjoint representation 
of SU C (3) algebra. In eq.(||), X q takes value | jj 3 - in the standard model, but here its value is 
irrelevant for our purpose. It should be mentioned that here the quark inserters are constructed 
by borrowing the fermion-Higgs- vertex of Yukawa type from the standard model, this is in analogy 
with as occurs in QED. The issue has been discussed in detail in 

For a given divergent 1PI amplitude T^ n f' Us \pi, ■ ■ ■ ,p„, f ; ki, ■ ■ ■ , k ng ) at the one loop order 
with nf external fermion lines and n g external photon lines, we consider a set of 1PI amplitudes 
p(«/,n 9 )^ 1) . . . p . ki, - ■ ■ , k ng ; q) which correspond to the graphs with, if the loop contained in 
the graph purely consists of fermion lines, all possible 2q insertions of the fermion inserter in the 



4 



internal fermion lines, or in other cases, all possible q insertions of the corresponding inserter in 
the internal boson (ghost) lines in the original graph. The divergent degree therefore becomes: 

S = 4-I f -2I g - 2q. 

If q is large enough, T^ n f' n ^(pi, ■ ■ ■ ,p nf ; k±, ■ ■ ■ , k ng ; q) are convergent and the original divergent 
function is the case of q = 0. Thus we reach a relation between the given divergent 1PI function 
and a set of convergent 1PI functions at the one loop order. In fact, the function of inserting the 
inserter (s) into internal lines is simply to raise the power of the propagator of the lines and to 
decrease the degree of divergence of given graph. 

In order to regularize the given divergent function with the help of this relation, we need to 
deal with those convergent functions on an equal footing and pay attention to their differences due 
to the insertions. To this end, we introduce a new function: 

r( n f' n ^{p 1 ,---, Pnf ;k 1 ,---,k ng ; q; fj,) 

(5) 

= (- i/U )^(-i\)- 2 ^ £r( n '> n »>(pi, . . . kl , ■ ■ ■ , k ng ; q) 

where // is an arbitrary reference mass parameter, the summation is taken over the entire set of 
such N q inserted functions, and the factor (— i\)~ 2q introduced here, in which A stands for \ q 
for fermion loop and for g c for other cases, is to cancel the ones coming from the inserters. It is 
clear that this function is the arithmetical average of those convergent functions and has the same 
dimension in mass, the same order in coupling constant with the original divergent 1PI function. 
Then we evaluate it and analytically continue q from the integer to the complex number. Finally, 
the original 1PI function is recovered as its q — > limiting case: 

r( n " n »)(pi, ••••/>,.,:/.•:••••. k ng ) = lim r( n /> n »)(pi, • • ■ ,p nf ;ki, • • • , k ng ; q;»), (6) 

and the original infinity appears as pole in q. Similarly, this procedure should work for the cases 
at the higher loop orders in principle. 

The divergent 1PI graphs at the one loop order in QCD are as follows: the gluon self-energy 
W^ v {k), the quark self-energy the ghost self-energy Tl ab (p), renormalized by Z3, Z% and Z3, 

the three-gluon vertex ^ a ^y\{k\, fo), the four-gluon vertex r^Ar(^ 1 ' ^ 3 )' ^ ne q uar k-gluon vertex 
r"(p',p), and the ghost-gluon vertex i"^ c (p',p), with the renormalization constant Z\, Z4, Zf, Z\. 
In addition, there is a mass shift for the quark, which we shall ignore. All corresponding graphs are 
listed in figures 1-7. As numerious diagrams are concerned, evaluating them one by one in detail 
would be much lengthy and unnecessary. In the next section, we will evaluate in detail the gluon 
self-energy IT^, (k) as a typical example to show the main step of the approach, paying special 
attentions to the gauge inva 
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riance of the function. Then in the subsequent section, we will directly give all the results 
corresponding to other involved diagrams to verify the Slavnov- Taylor identities at one loop order. 



3 Regular izat ion and Evaluation of the Gluon Self-Energy Il^(fc) 

Before the detailed evaluations are presented, we should first refer to a special problem which 
arises in any massless theories, i.e., the genuine infrared divergence in these theories. In the 
regularization schemes, this problem usually appears as the lack of consistent definitions of the 
regularized Feynman integrals for the ones which are both ultraviolet and infrared divergent. For 
instance, in the dimensional regularization scheme, let's consider the massive integral 

' — 7 = I(m, u,n), ( m 2 + ) (7) 



2\n—uj 



(2vr) 2w (I 2 + m 2 ) n (47r) w r(n)(m 

which converges for ui complex; the parameter n is arbitrary but fixed. We note first of all that 
the limit lim I(m, uj, n) may or may not exist, depending on the relative magnitudes of n and 

m 2 ^0 

uj. But even if it did exist, another problem could arise as we approach four-space (provided the 
original integral is infrared divergent to begin with), because in general 

lim [ lim I(m, uj, n)] ^ lim [lim /(m, uj, n)], 

m 2 -*0 ?n 2 ^0 

f d 2uj l 1 

so that the massless integral lim / - — ttt - ,, n , can not be derived unambiguously from the massive 

& u->2j (2vr) 2w (Z 2 )" & J 

integral (|^). Furthermore, the trick of inserting a finite mass into the integral and then allowing it 
to approach to zero at the end of the calculation is, in general, not a satisfactory prescription yet, 
because it spoils the gauge symmetry in the original theory, provided such a symmetry existed in 
the first place. To avoid this difficulty, 't Hooft and Veltman naively comjectured that 



d 2uJ l 1 
^2 J (2tt) 2 " (Z 5 ) 



lim o / /o-\2^ TmTT = ' foruj,n complex. (8) 



It has been shown that no inconsistencies occur, e.g., in the Slavnov- Taylor identities [10], due to 
the acceptance of the above conjecture. 

In our application of the present approach to QCD, as we will see, the same problem occurs, 
e.g., in calculating the tadpole diagram Fig. Id of the gluon self-energy. To solve this problem, we 
employ a conjecture analogous to 't Hooft and Veltman's: 

jW (2tt) 4 [{k-iy]M{p)Bq+n = °> f° r ?' n ' com P lex , A>0, B>0, A + B = l . (9) 
Likewise, we will see that no inconsistencies occur due to the acceptance of the eq.(S). 
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Now we turn to evaluate in detail the gluon self-energy Uf u (k) shown in Fig.l. The dia grams 
contributing to IT™, (k) are four in number, namely, the gluon loop diagram, the ghost loop dia- 
gram, the quark loop diagrams, and the gluon tadpole. The integral expressions of the regularized 
diagrams in the momentum space are given in the appendix ( For simplicity, we take the Feynman 
gauge £ = 1. ). 



First, we consider the gluon loop contribution. From ( 19a ), a little bit of algebra yields 

= -^[-6^(8)^]^/!, 

with 

T _ J_y- f d4 P ~ 2 P 2 9^u ~ (5k 2 - 2p ■ k)g^ u - Wp^p u + 2k^k u + 5p^k u + 5k^p u 
1 " N q f^J (2vr) 4 (p2)i+i[(p_A;)2]9-i+l 

In the present case, N q = q + 1. Note that because of eq.(^), the contribution of the first term in 
the numerator of the above equation actually vanishes, so it can be neglected. Using the Feynman 
parameterization, we get 



g + i^r(i + i)r(g-i + i)y 

d 4 p (2x - 5)k 2 9llu - \Qp^p u - (10x 2 - lOx - 2)k„k 



dxx q -Hl - xY 

rw _ 9M- h 

(10) 



(2tt) 4 [p2 + x{\ - x)k 2 ]i+ 2 

where we have made a momentun shift: p — > p — kx, and the linear terms in p in the numerator 
have been dropped since they do not contribute to the integral. Now the integration over p can be 
performed by using the following formulas: 

d 4 P (p 2 f i T(2 + P)T(A-2-P), 22+ p_ A 



(2tt) 4 {p 2 + M 2 ) A (At:) 2 T(A) 



(M 2 ) 2 +V- A , (11a) 



d 4 p (p 2 )%p„ i 1 T(3 + P)T(A-3-P) 2 3+0 _ A 

(2tt) 4 (p 2 + M 2 ) A ~ (4vr) 2 4 5 ^ r(A) [ ' ' { ' 



f d 4 p (p^fpftPvPpPa i 1 , . 

J (27r) 4 (p 2 + M 2 ) A ~ (47r) 2 2V 9>lu9pC7 9 ^P 9va + 3iur9vp) 



r(4 + p)r(A-4-p), 2 4+/3 _ A 



T(A) 

From @, (gig), and ( |TTb|) we get: 



(^g(g + l)^(^ ^ ^-^)^ 1 ( ) ,-^ + l)/o dxx ^ (1 xr 

(g - l)[(2x - 5)fc 2 ^ - (10x 2 - lOx - 2)k^k v ] - hx(l - x)k 2 g^ 



x 



[x(l - x)k 2 ]i 

* r(g-i) f 1 j [<l( 2x ~ 5 ) + ( 5x2 ~ 7x + 5)]fc 2 g^ - (g - l)(Wx 2 - Wx - 2)k^k u 
W7T^+2lJ ax [x(l-x)k 2 ]<i 
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where in the last step the summation over i has been performed with the help of the binomial 
theorem. This expression also makes sense when we make an analytical continuation of q from 
integer to complex number. When q — > 0, using the expansion 

Ifcil = _I _ q+ (q% [x(l - X )k 2 ]i = 1 + qHx(l - x)k 2 ] + 0(q 2 ) , 
T(q + 2) q 

and making a rescaling of the parameter [i 2 — > — 6C2(8)/i 2 , we get: 

n ab (k A _ ^C 2 (8) b 19 2 11 1 
n (yl)^l fc ^'MJ - K 9 ^ v y v v > q 

-C^g^u ~ y In(^) + (^V^ - yM^)] • ( 12 ) 

Clearly 11^^ (/c; /i) does not conserve current, this is due to our choice to use a covariant gauge 
(rather than, for example, an axial gauge). For the sake of this choice, we have to introduce 
spurious gluon polarization states. These spurious states must be removed by taking the ghost 
loop contribution into account. We could have computed with an axial or "ghost free" gauge, but 
it is usually much easier to use the simple Feynman gauge and add in the ghost contribution. 



To compute ghost loop contribution ( tL9b| ) , we use the same Feynman parameterization arriving 

at 

= - 5c 2 C 2 (8)[C 2 (8)^ 2 ]^/ 2 , 

with 

j j" 1 dx f d A p p^pv - x(l - x)k ll k p 



o 



(2tt) 4 [p 2 + x(l - x)k 2 ]i+ 2 



where momentun shift: p — > p — kx has been made, and the linear terms in p in the numerator have 
been dropped. After performing the integration over p and taking the limit q — > subsequently, 
we obtain 

n^C*;?;/*) = ^^^K^V + \ k M~ q 

1 1 /fc 2 \ 5 1 

-{^k 2 g^ v - -k^k v ) ln(-2 J + i^g^ + -k^k v )} . (13) 



Again, we find that IL nf' (k; q; /S) also does not conserve current. However, the non-conserving 



term in (|l3|) exactly cancels that in (|l2|), and makes n^^(A;; q; fi) +U^ B ^ u (k; q; fi) gauge invariant, 
namely, 

^tVjJU*;?;/*) + n (B j^(fc; g; n)] = . 



s 



Next, the quark loop contribution H^ c ^ u (k; q; fi) need not be recalculated since it is just 
H^ D (k; q; //), which can be found in |], |?]], multiplied by a color factor NfTr (T a T b ): 

n(oj*(*;9;M) = N f Tr(T a T b )U% ED (k;q^) = ^f/^ 5 a \k,k u - k 2 g, u ) 



(4vr) 5 

r 1 2 f 1 , , Xl A: 2 x(l -x)-m? 

x — + - - 2 / cbx 1 - x In i ^ + •••]. (14 

3g 3 Jo A* 



Obviously, H^ b u (k; q; fj.) itself is gauge invariant. 

Finally, as in dimensional regularization scheme, from eq.(^) we know that the contribution of 
the gluon tadpole diagram vanishes: 

U {D ^ u (k;q^)=0 . (15) 



From (^), (|l3[), (14), and (|l5[) we obtain the final expression of the regularized gluon self-energy 

liV ( 



nj*(fc;g;/x): 



n^(fc;g;M) = n {A $,(k; q; ax) + n (B) <£,(fc; g;/x) + U {G ^(k; q; ax) 

= TCT^V " W{[|C 2 (8) - |iV/C 2 (3)]i + •••}, (16) 

which is the same as that derived by means of other regularization methods and satisfies the gauge 
invariance condition k^H ab u {k] q; /x) = 0. 



4 Renormalization Constants and Slavnov- Taylor Identities at one 
Loop Order 

So far we have calculated the gluon self-energy IIj*(jfe) in detail by means of the intrinsic 
regularization method at one loop level. As we have expected, the result turns out to be gauge 
invariant. However, this is not enough for us to say that the method preserves the gauge invariance 
of QCD. As a complement, we should further show that the Slavnov- Taylor identities between the 
renormalization constants hold, namely, 

— = = £i = ?± (17) 
Z% Z F Z3 Zi 

It is these identities that guarantee that the renormalized theory possess the same gauge theory 
structure as the original one. Moreover, they are essential for proving the renormalizability and 
unitarity of the theory. 
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To verify these identities, all divergent 1PI graphs, not only U^ u (k), but also others, must be 
taken into account. After lengthy and tedious calculations, we get 

Zx = l + ^ |C7 2 (8) - ~JV/C2(3)]i (18a) 

Z, = l + [|cb(8) - ^JV/C 2 (3)]i (18b) 

^4 = 1 + ^ [™C 2 (8) - |jV/<7 2 (3)]i (18c) 

= 1 - (Jja[Cfe(8) + §Cb(3)]| (18d) 

^^"(^^ (186) 

^ = 1-A^- (18f) 



(4vr) 2 2 9 

g c 2 g 2 (8) l 
(4vr) 2 2 q 



z 3 = 1 + A^- (Ug) 



It is easy to see that the Slavnov-Taylor identities (|17]) indeed holds. For a close observation, 

we note that although these expressions are calculated in a specific gauge therefore their gauge 

dependence are not explicit, they are in fact all gauge dependent ( For instance, in the axial gauge, 

the effective Lagrantian has no ghosts and has the same structure as in QED, resulting in the 

identity Zf = Z% , which is patently untrue in the Feynman gauge.). We also remark that the 

fermion contribution to the vector boson quartic self-interaction must diverge if the Slavnov-Taylor 

i 

identities (|l7| ) are to hold because we explicitly see that the ratio Zf /Z% Z 3 2 , by which the bare 
coupling constant is related to the renormalized one, contains a fermion contribution. On the 
contrary, the corresponding box diagram of QED is finite. The reason is that in QED the box 
diagram's divergence vanishes only upon symmetrization of the external photon lines denoted only 
by their vector in 

dices, while in QCD the symmetrization of the external lines can be performed in two ways: 
by symmetrizing on both vector and group indices, which as in QED gives no divergence, or by 
antisymmetrizing on both vector and group indices. It is this new contribution which deverges. 
The same reason can be applied to the fermion contribution to the triple gauge vertex. 
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5 Concluding Remarks 



We have shown the main steps and results for the regularization of the divergent 1PI functions 
at the one loop order in QCD by means of the inserter proposal for the intrinsic regularization 
method. It turns out to be satisfactory: The gauge invariance is preserved manifestly and the 
results are the same as those derived by means of other regularization methods. It is natural to 
expect that this proposal should be available to the cases at higher loop orders in principle. 

The renormalization of the QCD under consideration in this scheme should be the same as in 
usual approaches. Namely, we may subtract the divergent part of the n-point functions at each 
loop order by adding the relevant counterterms to the action. The renormalized n-point functions 
are then evaluated from the renormalized action. In the limiting case, we get the finite results for 
all correlation functions. 

It should be mentioned that the method presented here is somewhat analogous to the analytic 
regularization method developed by Speer jy], 12|. However, the two methods are in fact different: 
As is well known, the analytic regularization method violates unitarity, this is due to the fact that 
it actually continuing the power of propagators in an arbitrary way. While in our method this not 
the case. Here we trie to find out a procedure of regularization from some physical principles. That 
is, giving a ultraviolet divergent process, one can always find a set of convergent function obtainable 
from existing Feymann rules and for the limiting case it turns to be the original ultra-divergent 
one. Or from another point of view, given a ultra-divergent function, one can always extract a 
set of convergent functions from a certain physical process, which is gauge invariant. After taking 
the limitation of the convergent functions, the divergent function is naturally regularized. There is 
nothing changed, the action, the Feynman rules, the spacetime dimensions etc. are all the same as 
that in the given QFT. From this viewpoint, one should have no doubt of gauge invariance and the 
unitarity of the method since the sum of the convergent functions comes from a certain physical 
process. 

Application of our approach to gauge theories containing spontaneous symmetry breaking such 
as the standard model should be straightforward. Also, it will be much helpful to apply our 
approach to some other cases, such as anomalies, SUSY theories etc., since in these cases the 
symmetries and topological properties are sensitive to the spacetime dimensions and the number 
of fermionic degrees of freedom etc., thus we are unable to tackle them consistently by means of 
the hitherto well-known regularization methods such as dimensional regularization method. It is 
reasonable to expect that the approach presented here should be able to get rid of those problems. 
We will investigate these issues in detail elsewhere. 
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Appendix: Integral Expressions of the Divergent 1PI Graphs at 
One Loop Level in QCD 



There are a number of divergent iPI graphs at one loop level in QCD, which contribute to the 
gluon self-energy U°^(k), the quark self-energy the ghost self-energy Tl ab (p), the three-gluon 

vertex T ab £ x {ki, k 2 ), the four-gluon vertex r^ c ^ r (/c 1 , fc 2 , k%), the quark-gluon vertex T^(p',p), and 
the ghost-gluon vertex T^ c (p',p) respectively. In what follows we present the integral expressions 
of the regularized diagrams in the momentum space (in Feynman gauge £ = 1). 

1. The integral expressions of the regularized diagrams contributing to II?* (k): 

1 2-2a . 2ar <<■ 2^ /oMom /T,ori6\ ^ f ^P 



n {A i(k;q;v) = ^ 2 ^ 2 "[-^g 2 c C 2 (8WTr(F a F b ) — £ / 
xG wiai (k, -p,p- k)G up2<72 (-k,p,k - p) 

/ — i \q—i+l 



(2vr) 4 



n (B t(fc;^) = <7rV 9 Hff c 2 ^2(8^Tr(F«F b ) ] l^| 

. p 2 / V (p — k) 



d A p 

W) 4 



1 i=0' 



n (c j*(fc;g;M) = ^ 29 HA s )^ / rr(T«T»)i-f;/ 



d 4 p 

(2^ 



/ 1 \ i+l / 1 \ 2o— i+1 



where Nf is the number of flavors, T a are generators of SU C (3) in fundamental representation, 
jabc j eno ^ e the structure constants of S77 c (3), (F a ) bc = —if abc are generators of SI/ C (3) in adjoint 
representation, and 

G fll/ \{pi,P2,P3) = {Pl-P2)x9^u + (P2-P3)^9uX + {P3-Pl)u9Xfi 
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comes from the three-gluon vertex. 



2. The integral expressions of the regularized diagrams contributing to rj^(fci, k 2 ): 

1 1-i , d A p 



1 i=0j=0 

x g p ^g^g p ^G vp2(J2 (k 2 , -k 2 - p,p)G Xp3r73 (-k 1 - k 2 , h -p, k 2 + p) 



r ( B&(*i, *2;</;m) = -^V'[-% c 2 c 2 (8)]'itj(|/' Sl11 

G ppiai {ki, -p,p - ki)[f bcr f t2S2r \g ua2 g Xp2 - g vp2 g\a 2 ) 



+f 2r f S2Cr (g» P2 g\* 2 - g v \g P2 o 2 ) + f 2r f 2r (gu\g P2 * 2 - g V a 2 g\ P2 )\ 

x**'^"v- 9 -(^) i+i ((^r i+i 

+ {(//, a, k x ) <-> (i/, 6, fc 2 )} + {0, a, fei) <-> (A, c, fc 3 )} , 

r ( c)t C A(fci, ^; </;/*) = - % 3 - 2 V 29 [-% 2 C 2 (8)] 9 Tr(F a F b F c )— \ T^TI 

9 i=0 j=0 ^ ' 

+ {(i/, 6, fe 2 ) <-> (A,c, fc 3 )} , 



ly Q i=o j=o J \ Z7T > 
x Tr [ 7 f * V +1 7A f * y' +1 7 i * \ 2g— i— j+1 

j^— jk\—m' V Jk 2 — m> v \fe — m> 
+ {(v,b,k 2 ) <-> (A,c,/c 3 )} , 

3. The integral expressions of the regularized diagrams contributing to r^ c A i r (/c 1 , fc 2 > fc 3 

rrf^i, fc 2 , k 3 ;q;n) = 5c 4 " 2 V 2<? [-6^ c 2 C 2 (8)]^r(F a F fe F c F d )^ £ £ *£f 
d 4 p 

x G rpi(T4 (-ki - k 2 - k 3 , ki - p, k 2 + k 3 + p) 

X G Ap?03 (A;3, -*2 ~k 3 -p,k 2 + p )gP^ g P2a 3g Psa 4gPi a 1 

x (5)' 



X 



/" d 4 p 



N 

ly 1 i=0 j=0 Z=0 



■ p 2 ) y(p+k 2 ) 2 ) y(p + k 2 + k 3 ) 2 ) V(p-/ci) 2 



+ {{v, b, k 2 ) -> (A, c, fc 3 ), (A, c, fc 3 ) -> (r, d, fc 4 ), (r, d, fc 4 ) -> (i/, 6, fc 2 )} 
+ {(i/,6,/c 2 ) ->■ (r,d, fc 4 )), (r,d,A; 4 ) -»■ (A,c, fe 3 ), (A, c, fc 3 ) — > (i/, 6, fc 2 )} 
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r ( B)#ft(*i. fa> fa;<z;/x) = - ^V«[-% c 2 c 2 (8)]^£ f 



(2^ 



P ct\9t P \ 9 PP \9ra\ to\ 9 P r9pi<Ji, 

+ f aSiri f dtiri (9pr9 Pl a 1 ~ 9^ l9 r Pl )][f bCr2 f t2S2r2 (9 V a 2 9X P2 ~ 9u P2 9Xa 2 ) 

+f^ n f 9 ^{9u P2 9x n -g v x9p a a 2 )+f^ 2n f^ n {9uX9p a a 2 -9u n 9xp a )] 
+ {(t, d, fa) <-> (A, c, fc 3 )} + {(r, d, fa) <-> (i/, 6, fa)} , 



(21b) 



JV 9j=0j=0 J 



d 4 p 

W) 



1 



x <y fli * 2 <y fl2 * 3 «5 S3ti [/ ad 7* lsir (5^ 1 9r P1 

l(Tl 9 P Ui9t Pi )\ 

x f bS2t2 f CS3t3 G up2a2 (k 2 , -k 2 -p,p)G Xpaa3 (k 3 , -k 2 -h-p, k 2 +p) ( 2 lc) 

9 9 9 V' {( p + k 2 f) {( p + k 2 + k 3 f) 

+ {(r, d, fc 4 ) <-> (A, c, fc 3 )} + {(r, d, fc 4 ) <-> (v, b, k 2 )} 
+ {(/x, a, fei) <-> (A, c, fc 3 )} + {(/x, a, fei) <-> (i/, 6, fc 2 )} 
+ {(//, a,ki) <-> {y, b, k 2 ), (r, d, fc 4 ) <-> (A, c, fc 3 )} , 



r (D )^(^i, fa, k 3 ;q^) = -^ V 2 g[ _ 6 ^2 C2(8)]gTr(F a F6FCF(i) _^^^^ J 



i=0 j=0 1=0 

jj-riP P {p + k 2 ) v {p + k 2 + fa) A (p - fei)r(— 



X 



(2vr) 4 ^ ""^ ' ±y ' V 

z i \ i+i / i \ <;-■> i i+i (21d) 



( i y + v \ y 

Wo + JfeoY 2 / \(n + ko + k?.) 2 J 



(p + k 2 ) 2 J V( p + k 2 + k 3 ) 2 J ^ip-h) 2 - 
+ {{v, b, k 2 ) -> (A, c, fc 3 ), (A, c, fe 3 ) -> (r, d, fc 4 ), (r, d, fc 4 ) -> (i/, 6, fc 2 )} 
+ {(i/, 6, fc 2 ) -> (r,d, fc 4 )), (r,d, fc 4 ) -> (A,c,& 3 ), (A, c, fc 3 ) -> (i/, 6, k 2 )} 
+ {(i/, 6, fc 2 ) <-> (A, c, fc 3 )} + {(i/, 6, fc 2 ) <-> (r, d, fc 4 )} + {(A, c, fe 3 ) <-> (r, d, fc 4 )}, 

r (B) 5? T (fci, fa, fa;<^) = -5 c V 29 (-^) 2 «Tr(T a T 6 T c T d )^X] E ' 

iV 9 i=0 j=0 Z=0 

x 7A ( * ) ,l -,(-^) 2 "' ' / "'; Pie) 



+ {(^ &, fa) (A, c, fe 3 ), (A, c, fe 3 ) -> (r, d, fc 4 ), (r, d, fc 4 ) -> (i/, b, k 2 )} 

+ {(u,b,k 2 ) -> (r,d, fc 4 )), (r,d, fc 4 ) -> (A,c,& 3 ), (A, c, fc 3 ) — > (i/, 6, fc 2 )} 

+ {(i/, 6, fc 2 ) <-> (A, c, fc 3 )} + {(i/, b, k 2 ) <-> (r, d, fc 4 )} + {(A, c, fe 3 ) <-> (r, d, /c 4 )} , 
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4. The integral expressions of the regularized diagram contributing to U ab (p): 

n ab (p;q;ri = -s c 2 - 2 V 2? H^ 2 (8)]'Tr(F«F^^. k (iy + \jk^) ' (22) 



5. The integral expressions of the regularized diagrams contributing to T abc {p' ,p): 

r^ip'^q-ti) = l gl-^^[- l glC 2 {S)fTr{¥ a Y b Y c )J-^k^p'P 

x Gfipcr (ki-k 2 ,p-k 1 M-p)(^) q+ \ {k _ % pr)2 )( {k ^ p)2 ) , (23a) 



/d A k 
T^u(k + p'Uk+ P )-p' 



6. The integral expressions of the regularized diagram contributing to T,(p): 
E(p; g;A1 ) = T a T a S^ D (p;g ;/ u) 



<i 4 /c z / — Z 



- -»^ 2 > 2, l- 6 ^^< 8 'l ,T ° T V(^ h, 'T^> ) (7^w)' +1 . < 24 > 



7. The integral expressions of the regularized diagrams contributing to T a ^(p' , p)\ 

xf ast G^(p' -p,p-k,k-p') 



(27T) 



4 



r ( B) a >',p;?;A*) = T 6 T a T fe r^V,p;<^) 

= -i^ 2 v» 2 *[-Gi£C 2 (8)]*T b T a T b J ' /:/ 



(2vr) 4 
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Figure Captions 

Figure 1 The one loop Feynman diagrams which are needed for calculation of the renormalization 
constant Z3. 

Figure 2 The one loop Feynman diagrams which are needed for calculation of the renormalization 
constant Z\. 

Figure 3 The one loop Feynman diagrams which are needed for calculation of the renormalization 
constant Z±. 

Figure 4 The one loop Feynman diagram which is needed for calculation of the renormalization 
constant Z%. 

Figure 5 The one loop Feynman diagrams which are needed for calculation of the renormalization 
constant Z\. 

Figure 6 The one loop Feynman diagram which is needed for calculation of the renormalization 
constant Z[ . 

Figure 7 The one loop Feynman diagrams which are needed for calculation of the renormalization 
constant Zf . 
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v, b, k 




Fig. la The gluon loop diagram. 
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\ u,b,k 



Fig. lb The ghost loop diagram. 



v, b, k 




Fig.lc The quark loop diagram. 




Fig. Id The gluon tadpole. 
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Fig. 3c 



Fig. 3d 
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